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ABSTRACT 

In this paper, we show that the supernova remnant (SNR) masses cumulated 
from core-collapse supernovae along the star formation history of two powerful 
z=3.8 radio galaxies 4C41.17 and TN J2007-1316 reach up to > 10® Mq, com¬ 
parable with supermassive black hole (SMBH) masses measured from the SDSS 
sample at similar redshifts. The SNR mass is measured from the already ex¬ 
ploded supernova mass after subtraction of ejecta while the mass of still luminous 
stars fits at best the observed spectral energy distribution (SED), continuously 
extended to the optical-Spitzer-Herschel-submm domains, with the help of the 
galaxy evolution model Pegase.3. For the recent and old stellar populations, SNR 
masses vary on 10®“^® Mq and the SNR-to-star mass ratio between 1 and 0.1 % 
is comparable to the observed low-z SMBH-to-star mass ratio. For the template 
radio galaxy 4C41.17, SNR and stellar population masses estimated from large 
aperture (>4arcsec=30kpc) observations are compatible, within one mass order, 
with the total mass of multiple optical HST ( 700pc) structures, associated with 
VLA radio emissions, both at O.larcsec. Probing the SNR accretion by central 
black holes is a simple explanation for SMBH growth, requiring physics on star 
formation, stellar and galaxy dynamics with consequences on various processes 
(quenching, mergers, negative feedback) and a key to the relation bulge-SMBH. 
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Subject headings: galaxies: evolution—galaxies: high-redshift—galaxies: individ¬ 
ual (4C41.17, TNJ2007-1316)—(galaxies:) quasars: supermassive black holes— 
(cosmology:) early universe 


Introduction 


The scaling relations found up to z = 5 flMaeorrian et ah Il998l: iFerrarese &: Merritt 


20001) between the supermassive black hole mass Msmbh and the properties of their spheroidal 
host such as the stellar velocity dispersion a, the 17-band luminosity L an d the bulge stellar 


mass Mbuige have been recently conhrmed at the center of nearby galaxies flMcConnell et ah 


20131 ). However, the origin of such a tight relationship between SMBH and galaxy star for¬ 
mation history is still actively debated. The main question is to justify the huge difference 
of size scales between the galaxy core (~ Ipc) embedding the supermassive black hole, and 
the galaxy and its environme nt (from 10 to 100 kpc). At the core scale, the central black 
hole grows by fluid accret ion (IBlandford fc Begelman Il2004r) . and already reaches 10®^^° Mq 


at the highest redshifts fiVestergaard et ah 1120081) . From the gas point of view, because 


disk-outflows are estimated unsufficient to justify the SMBH growth, AGN feedback is si¬ 
multaneously proposed to expl ain the Msmbh-ct relation and to quench the star formation 
activity in early type galaxies flGabor fc Bournaud II2014J) . However, even if a large number 
of feedback models are supported by many surveys (X-rays, radio VLA, PEP, GOODS- 
Herschel, GANDELS-Herschel), difficulties remain to explain any accretion process from 
variable AGNs, which have never been observed to be turned on or off by star formation. 

From the star point of view, detailed analyses of stel l ar processes n ear the massive 
black hole in the Galactic Genter system (lAlexander l2005l: iGenzel 120141) study the mass 
segregation driving stellar dense matter (neutron stars and ste l lar black holes) towards th e 
central black hole by dynamical friction flGhandrasekhar Ill943l: iBinnev &: Tremaine 1120081) . 


Adopting the universality of the IMF, t he SNR mass of old ste llar populations is waited for 
a few percent of the total galaxy mass flMagorrian et ah 1119981) . However at high redshifts. 


stellar population properties (mass, size, density, velocity dispersion, drag time-scale) depend 
on cosmology and galaxy evolution: at z=3.8, the lookback time is more than 80% of the 
cosmic time. So that if SNRs are eternal and have time to segregate, the SNR-black hole 
relation by dynamical friction depends on galaxy evolution parameters, in particular the star 
formation rate (SFR) by galaxy types, and evolving star populations need to be followed 
with an evolutionary code as the code Pegase. 

The AGN-starburst relation is often based on star formation rate estimates derived 


mostly from mid- and far-infrared luminosities flKennicutt 1119981) . However, the SFR cal- 
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ibration suffers from uncertainties in the disentangling of the emissions by the AGN torus 
and by cold grains, in the radiative transfer modeling of grain absorption/emission through 
various geometries and in the ultraviolet to optical contribution of young star photospheres. 
Moreover, a time-delay of the dust emission peak is induced by the lifetime duration of super¬ 
novae, the ejecta of which are sources of the interstellar metal enrichment and dust grains. As 
a consequence, most SFR estimates derived from the observed far-infrared emissions are due 
to less massive populations, still luminous after the supernova explosions, implying underesti¬ 
mation to be corrected with evolution models. The SFR calibration needs coherence between 
various tracers by taking into account stellar passive evolution w ith the help of evolutionary 
spect rophotometric and chemical models, e.g. the Pegase.3 code fjFioc fc Rocca-Volmerange 


19971 ). recently improved with grain C and Si enrichment and radiative transfer computed 


by Monte-Carlo simulations (near submission). 


Powerful radio galaxies are preferentially selected for their detection at the highest red- 
shifts (z>5), for their signatures of starburst and AGN activities through continuum and 
emission lines and for the presence of a supermassive black hole. From the catalog Herge 


flDrouart et ah I2014J ) of 70 1 < z < 5.2 radio galaxies observed with the Herschel and 


Spitzer satellites, completed by UV, optical and submm observations, the two z=3.8 radio 
galaxies 4C 41.17 and TN J2007 — 1316 have been selected for their low AGN activity fa¬ 
voring the domin ant stellar origin of the far-infra red dust emission. Their continuous SEDs 
have been htted flRocca-Volmerange et ah 1 120131) according to t he usual evolutionary spec¬ 
tral synthesis procedure fjhe Borgne fc Rocca-Volmerange 2002). Fits are processed in the 
observer frame (and not in the rest-frame) using libraries of template spectra for starbursts 
and Hubble Sequence galaxies, corrected for universe expansion and galaxy evolution, with 
various input parameters (initial metallicity and mass function IMF). 


Hereafter the new point is to estimate the SNR masses predicted by star formation his¬ 
tories of elliptical galaxies supposed to host distant radiogalaxies, assuming the universality 
of the IMF. Derived from evolutionary models, calibrated on low resolution SEDs, SNR and 
star masses are globally compatible with the high resolution (0.1 arcsec) HST and VLA 
data. Galaxy evolution models follow at any galaxy age and type i) the already exploded 
supernova population; ii) the consequent metal-enrichment; iii) the SNR mass ; iv) the mass 
of residual luminous stars. In section 2, the SNR mass of a Single Stellar Population M|^^(t) 
and of Hubble Sequence galaxy scenarios M^p^(t) are derived at age t. Section 3 shows the 
best-fits of SEDs, stellar masses, metallicities and SNR masses for the two radio galaxies. 
Section 4 analyses the high resolution images of 4G41.17 in terms of SNR and resolved stellar 
masses and of their compatibility with radio and optical (700pc) structures. The last section 
is a discussion on the observed tight bulge-black hole mass relation as well as on model 
degeneracy and perspectives. 
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2. Evolution of cumulated masses of SNRs and stars 
2.1. Single stellar popnlation (SSP) 

In the SSP formalism, the number of stars formed per mass and time units is at time t: 


(PN{m,t) = 6{t) $(m) dm dt 


where 5(t) is the instantaneous star formation rate SFR{t) and <h(m) is the initial mass 
function (IMF). 

The SNR (black hole or neutron star) is produced during the terminal phase of the iron 
core-collapse of massive supernova explosions: SNIIe, SNIb and others, followed by a shock 
wave rebound ejecting the envelop of e nriched gas. At any age, the superno va number and 


their corresponding ejecta from yields fjPortinari et ah Ill998l: iMarieo Il200ll) . depending on 


the initial metallicity of the SSP, are followed with the help of the evolutionary code Pegase 
flFioc fc Rocca-Volmerange 1119971) (see also www2.iap.fr/pegase) in its last version Pegase.3 
(near submission). 

At the age t of the SSP, the supernova remnant SNR mass ttisnr is cumulated from the 
successive generations of stars exploding as supernovae between 0 and t: 


/IfSSP 


it) = 


r.min(t™,t) 


'ms NR i'mpsN i'mpsN {i')) 


I 4-/ _ -f-rnin 

^ ~^PSN 


dm 


PSN 


dt 


{t')dt' 


where msNnijnpsNii')) = ijnpsN —'mej){t') with mpsN the presupernova mass, m^j the mass 
fraction of metal-enriched gas ejecta, including winds during evolution, and tpsN its lifetime 
duration depending on initial mass exte nded from the min up to the max progenitor lifetime 
durations. The Kroupa standard IMF flKroupa II1997I) is limited between the extreme values 
Tfsnp=120 Mq and 0.09Mq with supernova precursors from I 2 OM 0 to SMq. The 

stellar mass Mgt®^(t) of the SSP at the age t corresponds to masses of still luminous stars, 
less massive than the past supernova progenitors and passively evolving according to the 
stellar evolution tracks. The respective evolutions of the current Mg.^®^(t) and of 
are opposite up to ~ 50 Myr, the lifetime duration of the lowest-mass supernovae. 


2.2. SNR and stellar populations derived from galaxy scenarios by types 

The star formation laws by galaxy type modeled with the code Pegase.3 (Fig. [1], top) are 
calibrated on local templates of the Hubble Sequence by htting only three gas-regulation pa- 
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Fig. 1.— Top: Adopted star formation laws for types Ell, SO, Sa, Sd and Irregula r (Im) 
galaxies as a function of z in arbitrary units (fbe Borgne fc Rocca-Volmerange 2002 b Bot¬ 
tom: Comparison of synthetic templates of Ell (red), SO & Spirals (turquoise), Im(dark blue), 
quenched (pink) galaxies with observed loc al SDSS galaxies (black points) in the color-color 


g-r/r-i diagram flTsalmantza et ah 
formation laws 


20121) ■ Green points correspond to the adopted star 
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rameters (infall time scale, outflow age, gas-dependent star formation rate as the Kennicutt- 
Schmidt law) and the choice of the IMF (see www.iap.fr/pegase.2). The comparison with 
the observed local SDSS data sample is presented in the color-color g-r/r-i diagram (Fig. 
m bottom). Excepted for the most extreme red and blue galaxies, star formation laws give 
synthetic templates compatible with observations. Laws by types mainly differ by their 
star formation time scales from ~ 1 Gyr for early-types up to ~ 10 Gyr for late spirals 
flRocca-Volmerange et ah 112004) . In particular, the elliptical galaxy scenario shows its max¬ 
imal efficiency at z ~ 8 for a formation redshift Zfor=^0- Because of the cosmic time-redshift 
relation t-z, changing Zfor in the interval [5, 30 or more] will not signihcantly modify the 
scenario parameters. For each galaxy type, the SNR mass M^^(t) is derived from the 
cumulated number of supernovae, after subtraction of ejecta. 


3. Masses and ages of the 2 ; = 3.8 radio galaxies 4C41.17 and TN J2007-1316 


T he SNR and star masses are deri ved from the best hts of the continuous UV to submm 
SEDs dRocca-Volmerange et ah 2013), covering Herschel, Spitzer, VLT and submm ground- 
based data. The choice of faint AGN radio galaxies allow to disentangle the emission of 
stellar populations from that of the AGN-torus (FiglJl green dashed line). The procedure 
is based on large libraries of SSPs and galaxy scenarios by types built at ages 0 to 20 Gyr, 
various IMFs (Kroupa and Salpeter), initial metallicities {Z = 0 to 0.02) and gas-and-metal 
column density mass/r^ varying as 0.1, 1, 10 x the local gas column density NHL 


For the two radio galaxies, the SSP model hts the optical and far-IR data (Fig. [2l blue 
line) at the same age of 30-35Myrs (the lifetime duration of a 9Mq supernova) with Kroupa 
IMF in a dense medium (10 x local value ~ lO^^cm”^). The old elliptical/SO model mainly 
hts the SED peak in the near-IR to Spitzer domain, at cosmic ages ~ 1 Gyr (respectively 
0.7 Gyr and 1.2 Gyr) with Kroupa IMF in a high-density medium (Fig. [21 orange line). 
In agreement with ages and masses of luminous stars from these best-hts, SNR masses are 
derived for the young and old past stellar populations of the two galaxies. 


Table [H provides the age, star mass Mgtar, SNR mass Msnr, the MsNR/Afstar ratio and 
metallicity Z of the young (evolved SSP) and old (Elliptical/SO) stellar components for the 
two radio galaxies, both best-htted with the respective values of the reduced = 0.96 for 
4G41.17 and y^ = 1.53 for TN J2007-1316. SNR masses vary from 3.4 to 19 x IO^Mq with 
ratios MsNR/Tfstar varying from 2.1 x 10“^ to 4.2 x 10“^. The striking result is that the 
supernova remnant mass Msnr is of a few 10® to 10 ^®Mq, quite comparable to supermassive 
black hole masses. SNR masses Msnr increase by less than a factor 2 when changing from 
the Kroupa IMF to the Salpeter IMF: the massive star slopes are <h(m) oc (Kroupa) 
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Fig. 2.— Observations are red crosses and global SED fits are black lines for the two galaxies. 
Left: 4C41.17, yonng starbnrst (bine line) with age = 30Myr, M^tar = 2.2 x 10^^ Mq, Z = 
510“^ pins an early type old scenario (orange line) with age = 0.7 Gyr, Mgtar = 2.0x10^^ Mq, 
Z = 5.7 X 10“^. Right: TN J2007-1316, yonng starbnrst (bine line) of age = 35Myrs, 
Mstar = 0.8 X IO^^Mq, Z < 10“^ and early type old scenario (orange line) of age = 1.2 Gyr, 
Mgtar = 9.4 X IO^^Mq, Z = 1.55 x 10“^. Black-body or Pier AGN models (dashed green 
line) are compatible with snperior limits (green arrows). Standard cosmology parameters 
(Gm = 0.3, Ga = 0.7, Ho = 72 km s~^ Mpc~^). 




















and oc m (Salpeter). Exotic IMFs, as top-heavy ones, are not considered as they are 
unable to predict the typical 1 fim peak of red population for evolved early-type galaxies. 


The supernova remnant masses Msnr ar e quite comparable to sup ermassive black hole 
masses Msmbh measured at similar redshifts flVestergaard et ah 11200811 from the SDSS DR3 
15 000 quasar sample with redshifts between 0.3 and 5.0. At z = 3.8, the medium value is 
10®'^ in excellent agreement with our estimates of SNR masses Msnr- 


To summarize, the evolved luminous stellar population fitting the SED of z=3.8 radio 
galaxies allows to estimate a population of already exploded supernova remnants, the mass 
MsNR(t) of which is of the order of supermassive black hole masses Msmbh observed at 
similar redshifts. The same result is found for the two distant radio galaxies, both hosted 
by massive ellipticals. 


4. SNRs and stellar populations at low and high angular resolution 


SNR masses Msnr of 10® Mq, comparable to the 10®'^ Mq of SMBH masses from 
the SDSS3 Quasar sample, also agree with galaxy masses derived f rom the VLT broad H a 


line emission from the nucleus of distant powerful radio galaxies fiNesvadba et ah Il201lh . 
Moreover from the best fits, the ratio MsNu/Mgtars is to 10“^ c omparable to the lo cal 


ratio of MsMBn/Mstars- Recently the discovery of the Himiko galaxy fjOuchi et ah II2013I) at 
z~7, with low dust and metal content, would reveal the preliminary phase of intense star 
formation through UV light and CII line, preliminary to the supernova explosions and to 
dust attenuation from their ejecta enrichment. The current star formation rate of the evolved 
component of the two radio galaxies is respectively 243 Mq yr“^ and 77 Mq yr“^ giving an 
SFR surface density of a few hundred of Mq yr“^ kpc“^ if the stars are distributed over a 


radius of about 750 pc around the centre at age of ~ IGyr flWalter et ah 1120091) . 


Galaxy /component 

Age 

Star Mass 

SNR Mass 

MsNu/Mstar 

Z 

4G41.17/ SSP 
TNJ2007-1316/ SSP 

30 Myr 
35 Myr 

2.2 X 10^^ 
8.1 X 10^® 

8.5 X 10® 
3.4 X 10® 

3.8 X 10"® 
4.2 X 10-2 

0.0005 

0.0001 

4G 41.17/ SO 
TNJ2007-1316/ Ell. 

0.7 Gyr 
1.2 Gyr 

2.0 X 10^^ 
9.4 X 10^1 

4.3 X 10® 
1.9 X 10^® 

2.1 X 10-® 
2.0 X 10-2 

0.0057 

0.0155 


Table 1: Age, star mass M^tar , SNR mass Msnr, MsNR/Mgt^r ratio and metallicity Z of 
the young (SSP) and early-type (elliptical/SO) components of the two 2 ; = 3.8 radio galaxies 
4C41.17 and TN J2007-1316. Mass unit is in Mq. 
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Our results depend on the adopted evolution scenario of elliptical galaxies which is ro¬ 
bust. Forming the bulk of stars within a 1 Gyr time-scale, this scenario predicts the reddest 
colors of ellipticals in the z=0 color-color diagram, traces the shar p limit of the z=0 to 4 K-z 


Hubb le diagram with the evolution of radio galaxy hosts flRocca-Volmerange et ah 


20041) on £g. 4, reproduces photometric redshifts, comparable to spectro scopic ones, up to _ 

redshifts z=4: and justify the faintest mu lti-wavelength deep galaxy counts (iLe Borgne fc Rocca-Volmerange 
2002; Fioc fc Rocca-Volmerange 1999 ). 


Another factor is the angular resolution of 7.23 kpc arcsec“^ at z=3.8 with the cos¬ 
mology Ga = 0.73, Gm = 0.27 and Hq = 71kms“^Mpc“^. At low spatial resolution and 
for the two radio galaxies, SNR masses a re derived from the continuous optical-IR-submm 
SEDs rebuilt at the aperture ~ 4 arcsec flRocca-Volmerange et al. 1120131) . At about a sim- 
i lar aperture ^ 3-4 arcse c , SMBH masses are measured through SDSS broad-band hlters 
(jVestergaard et al. 1120081) . SNR and SMBH masses are measured on the similar large scale 
of 20-30 kpc. At hi gher resolution, th e template radio galaxy 4C41.17 was observed with 
HST in the optical (IMilev et al. Ill992h and w ith VLA on multifrequency radio continuum, 
in particular at 8.3 GHz flGarilli et ah. 1119941) . At the same 0.1 arcsec resolution, multiple 
structures of ~ SOOpc size are identified. The remarkable similarity between the optical and 
radio morphologies of these structures is conhrmed: the evolved stellar components as HI to 
H4 are within the error bars, located on compact radio emission zones B1-B 2-B3 with the 
core N between B1 and B2, see hgure Ic and plate 55 of (IGarilli et ah. Ill994l ). The similar¬ 
ity of the individual structures implicitely link evolved stars and radi o jets. More recen tly 
near-infrared integral held spectroscopy with Gemini ALTAIR NIFS flSteinbring Il2014j) at 
the similar resolution of 0.1 arcsec conhrms the clumpy evolved structures and identihes a 
bow shock, possibly tracing recent star formation. Interpreted as active nuclei of galaxies 
from their hat radio emissions, the 0.1 arcsec resolved images witness the old evolved stars 
and black holes emitting radio jets, possibly fueled by SNRs. The rough approximation of a 
compac t black hole population from a rapid st arburst episode of 10®yr due to jet-triggered 

' TIE ' 


shocks fiBicknell et al. Il2000l: ISteinbring 1 120141) . spiraling towards the center of the galaxy 


throu gh the stellar structure supposed a singular isothermal sphere flBinnev &: Tremaine 


20081) . loosing its angular momentum by dynamical friction to the galaxy core, gives the 
following time scale : 

IO^Mq 


Ri 


19Gyr 


Tic,clump 


(ftT) 


a 


= 9.510-^Gyr 


Ink ^5kpc 200 km s ^ Msnr 


for Ti = 0.25/cpc, Msnr = 2 x IO^Mq of one indivi dual structure with MsN R/M,^tnr.‘<=2 x 10 


with the highest observed value a = 400fcms“^ flMcGonnell et al. 1120131) and a low impact 
parameter Goulomb logarithm Ink = ln{bmax/b90) = 5. This average value of ~ 0.01 Gyr 
respects the Hubble time. However more rehned computations are required because mass 
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and number of stellar black holes are depending on the starburst age and IMF. 


The high resolution observations confirm the 500-700 pc structures witness a interme¬ 
diate phase of evolution before building up of exponential disks, dense bulges and super- 
massive black holes. At z=3.8, the evolution time-scale of ~ IGyr is too short, a longer 
evolution time-scale of secular evolution is cl e arly needed w ith the help of numerical models 


flNoguchi. M. I 1999 I: lElmeereen et ah 120081: iGenzel 1201411 but respecting the ^=3.8 con 


straints. The SINFONI/VLT data on z ~ 2 star-forming galaxies f Genel et ah II2OO8 : Genzel 


2 OI 4 J ) suggest a scenario predicting the fraction of the central bulge mass increasing secularly 
over timescales of less then a few Gyr, even without rnaior mergers. Recent observations of 
dense cores in massives galaxies fjVan Dokknm et al. II2014J) out to z = 2.5 are compatible 
with these early-type evolution scenarios. 


5. Discussion 

The Herschel far-IR emission is not as a direct tracer of star formation rate as are ionized 
gas emission lines or the UV-light continuum: a time delay is required to enrich the ISM in 
O, G and Si through supernova ejecta and to produce the dust grain mass required to £t 
Herschel data. Typically, the maximum far-IR emission peaks at age ~ 1 Gyr in the elliptical 
scenario. 

After such a time delay, supernova explosions created a population of dense matter 
remnants (stellar black holes and neutron stars), the number and mass of which depend on 
star formation history (rate and IMF) and galaxy age while the still alived luminous stellar 
population is measurable through the multi-wavelength galaxy SED. 

The accuracy of SNR and bulge masses depends on the quality of the minimization 
procedure fitting SEDs, mainly constrained by the far-IR/submm luminosity coherent with 
the UV/optical attenuation factor while the mid-IR may reveal the old giant star population 
at z=3.8 (see Fig2. bottom: the Imicron peak of the red sequence in TN J2000-1316 SED). 
The degeneracy of modeling, higher (60%) for the old than the young one’s (20%) is however 
limited because the optical domain is strongly contrained by the attenuation by metals in 
coherence with the far-IR. Far-IR is essential but not sufficient for the young component 
and remains less constraining for the old components. The AGN model and the far-UV 
polarisation (3%) which could contribute to the degeneracy are faintly concerned in the 
analysis of these two radio galaxies. 

Finally, owing to the high spatial resolution, individual stellar structures or clumps 
associated to radio jets, are identifyied during a so called intermediate phase of galaxy 
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evolution. In each structure, the bulk of evolved stars is already formed from intense star 
formation rate within ~ IGyr, its SNR population would migrate towards the nucleus fueling 
the central black hole, source of the radio jet. The further phase, not yet seen at z = 
3.8, is the formation of dense massive bulges and SMBH of ellipticals by coalescence of 
these individual clumps with their massive black holes and stars without need of extra 
formation of dark matter. The implementation of this intermediate phase of early-type 


galaxy evolution in detailed numerical simulations o f bulge formation in galaxies fISales et ah 


2OI2I : lElmegreen et ah I2OO8I: IPubois et ah 


2014 1 or in clusters f Antonini 2014j) is then 


required. 


An important consequence of the migration of the totality of SNRs formed during the 
evolution scenario of early type galaxies is to grow SMBH as star masses contribute to dense 
massive bulges and likely dnst tori in radio galaxies. No more need of extra formation of 
dark matter by black hole mergers or qnenching star formation or negative AGN feedback, 
all are complex processes not significantly confirmed by observations. More modeling on 
stellar dynamics in agreement with star formation process in coherence with galaxy evoltion 
models is then reqnired 


Statistical results from SMBH mass functions compared to stellar mass fnnctions at 
high redshift, in particnlar from large snrveys as EUGLID will dehnitely conhrm the keys 
of the AGN-star formation relations. Moreover the transfer of SNRs and inspiraling process 
of snch a hnge nnmber of stellar black holes contribnting to the snpermassive black hole 
growth, as possible cases of coalescence of two nnclear galaxy black holes, will be among the 
most i mportant class of gravitation al wave sonrces, for the fntnre space-based interferometer 


eLISA f Amaro-Seoane et ah 11201211 . 


We thank the referee for the helpfnl remarks that have improved the mannscript. We also 
thank Michel Fioc for many developments in the code Pegase.3 and Jean-Philippe Beanlien 
for his well-advised comments. This work is based on observations made with the Herschel 
Telescope, an ESA space observatory with science instrnments provided by Enropean-led 
Principal Investigator consortia and with important participation from NASA, and with 
Spitzer Space Telescope which is operated by the Jet Propnlsion Laboratory, Galifornia 
Institnte of Technology nnder a contract with NASA. 
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